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ABSTRACT 

Context. This is the second paper of a series devoted to the Wide Field Nearby Galaxy-cluster Survey (WINGS). WINGS is a long term 
project which is gathering wide-field, multi-band imaging and spectroscopy of galaxies in a complete sample of 77 X-ray selected, 
nearby clusters (0.04 < z < 0.07) located far from the galactic plane Qb\ > 20"). The main goal of this project is to establish a local 
reference for evolutionary studies of galaxies and galaxy clusters. 

Aims. This paper presents the optical (B,V) photometric catalogs of the WINGS sample and describes the procedures followed to 
construct them. We have paid special care to correctly treat the large extended galaxies (which includes the brightest cluster galaxies) 
and the reduction of the influence of the bright halos of very bright stars. 

Methods. We have constructed photometric catalogs based on wide-field images in B and V bands using SExtractor. Photometry 
has been performed on images in which large galaxies and halos of bright stars were removed after modeling them with elliptical 
isophotes. 

Results. We publish deep optical photometric catalogs (90% complete at V ~ 21.7, which translates to ~ M v + 6 at mean redshift), 
giving positions, geometrical parameters, and several total and aperture magnitudes for all the objects detected. For each field we have 
produced three catalogs containing stars, galaxies and objects of "unknown" classification (~ 16%). From simulations we found that 
the uncertainty of our photometry is quite dependent of the light profile of the objects with stars having the most robust photometry 
and de Vaucouleurs profiles showing higher uncertainties and also an additional bias of 0.2"'. 

The star/galaxy classification of the bright objects (V < 20) was checked visually making negligible the fraction of misclassified 
objects. For fainter objects, we found that simulations do not provide reliable estimates of the possible misclassification and therefore 
we have compared our data with that from deep counts of galaxies and star counts from models of our Galaxy. Both sets turned out to 
be consistent with our data within ~ 5% (in the ratio galaxies/total) up to V ~ 24. 

Finally, we remark that the application of our special procedure to remove large halos improves the photometry of the large galaxies 
in our sample with respect to the use of blind automatic procedures and increases (~16%) the detection rate of objects projected onto 
them. 

Key words. Galaxies : Clusters : General - Catalogs 



1. Introduction 

Clusters of galaxies are privileged systems to study, among oth- 
ers, two basic problems. On one hand, they (or a substantial 
part of them at least) are the largest gravitationally bound struc- 
tures and as such, physical entities whose properties must be ex- 
plained by cosmological theories. On the other hand, clusters are 
made of hundreds or even thousands of galaxies in high density 
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* Based on observations taken at the Issac Newton Telescope (2.5m- 
INT) sited at Roque de los Muchachos (La Palma, Spain), and the 
MPG/ESO-2.2m Telescope sited at La Silla (Chile). 



regions whose properties can be studied and straightforwardly 
compared and their evolution analyzed. 

Clusters of galaxies were first detected as marked overden- 
sities in the projected number of galaxies. Even if today it has 
been recognized that galaxies represent only a small fraction of 
the total mass of the clusters, they still are a fundamental tool to 
study cluster properties since they are usually much more easily 
detected and measured than the X-ray emitting intracluster gas 
or the more evasive non-baryonic component. 

Clusters of galaxies also have been widely used to study 
the evolution of galaxies in dense environments. Some of the 
first clear evidence of evolution of the galaxies came from 
the study of the populat ions of galaxies in clusters made by 
iButcher & Oemlej dl978h . They found that the fraction of blue 
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galaxies was higher in clusters at z > 0.4 than in nearby clus- 
ters, which they interpreted as the result of the aging of spi- 
ral galaxies after losing their gas supply and, therefore, dimin- 
ishing their star formation rates. In parallel, iDresslerl dl980l) 
showed that the central, denser parts of clusters are mainly pop- 
ulated by early type galaxies. Since then, several works have 
discussed the idea that morphological c ontent in clus t ers of 
galaxies changes along the Hubble time dDressler et ail 1997 



Fasano et all 120001: Ivan Dokkum et all 120001: iLubin et alll2002 ; 



Postman et ail 120051) . 

Two are the main reasons that make clusters useful to study 
the evolution of galaxies. First, galaxies in clusters can be con- 
sidered as being at the same distance since usually this dis- 
tance is much greater than the linear dimensions of the clus- 
ters. Therefore, knowing the redshift of a small subsample of 
the cluster is enough to (statistically) know the distance (once 
background correction it is allowed for) of hundreds or even 
thousands of galaxies. Of course, this highly increases the statis- 
tics when analyzing the properties of the galaxies. Second, clus- 
ters of galaxies can be detected up to high redshifts using dif- 
ferent techniques that are sensitive to local enhancements of the 
galaxy density or through the X ray emission from the intraclus- 
ter medium. 

Using the first all sky survey (The Na tional Geogra phic 
Society-Palomar Observatory Sky Survey, lAbelll 1 1 9591) the 
first sys tematic catalogs of clusters of galaxies were con- 
structed (lAbelllll957T:IZwickv et al.Lll963l) . However, it was only 
in the 1970s when clusters of galaxies started to be used as 
laboratories to study the properties of the galaxies they con- 
tain and how their evolution is af fected by the environment 
dGunn & Gotdll972tlOemlerlll974 . From this point on, the in- 
terest in clusters of galaxies at always higher redshift continu- 
ously grew and nowadays, specially thanks to the Hubble Space 
Telescope, clusters at z > 1 hav e been studied in de tail (e.g. 
RDCS J 1252 2927 at z = 1.235; iPostman et all l2005h and us- 
ing differejittechnique^ > 2 
(e.g. lSteidel etaill200alKurk et alll2004l) . 

To produce a correct interpretation of the observations at dif- 
ferent redshifts a good knowledge of the cosmic variance of the 
properties of clusters and galaxies at each redshift is needed, 
to ensure that the changes observed with distanc e are statisti- 
cally s ignificant. Paradoxically, as we discussed in lFasano et alj 
d2006l hereafter Paper I), there is a relative lack of knowl- 
edge of the properties of the clusters of galaxies in the local 
Universe at z ~ and the comparison of higher redshift clus- 
ters is always done with single clusters such as Virgo, Coma 
or Fornax. Indeed, for more than 20 years the most complete 
study of galaxie s in nearby clusters of galaxies has been that by 
IDresslerl (1 1 980t> . Dressler's work was based on data extracted 
from photographic plates. Several programs have been set up 
to continue that work with modern tools and techniques. Each 
of them addresses a particular aspect or problem and, there- 
fore, we still lack a general local reference for evo lutionary stud- 
ies. Thus, the ESQ Nearby Abell Cluster Survey dKatgert et all 
119961: iBiviano et ail I1997L ENACS) includes spectroscopy of 
galaxies in a number of clusters but the imaging is not as 
deep as that of Dressler's work. The Las Campanas/ATT R ich 
Cluster Survey dO'Helv et all Il998fc iPimbblet et atl 1200 lb is 
deeper but includes only 20 clusters. I n rece nt years, large sky 
surveys as the 2dF ( De Prqpris et all 2002 ) and , particularly , 
the S PSS dGoto et all I2002t iBahcaU et all 120031: iMiller et all 
120051) have been sources for the compilation and analysis of 
large samples of clusters of galaxies. However, while the 2dF 
is only spectroscopic the SDSS, which covers a very large area, 



is not deep enough to study the faintest part of the luminosity 
function. More rec ently, the NOAO Fundamental Plane Survey 
dSmith et all 120041) started to fill this lack of data, however its 
main goal is not the study of the evolution of galaxies and clus- 
ters but the large scale velocity fields using the Fundamental 
PlaneQ 

In t his con text, the Wide-field Nearby Galaxy-clusters 
Survey dPaper J) has been presented with the specific goal of 
sampling the properties of clusters and galaxies in clusters in 
the local Universe. This means to establish both the average and 
the variance of the properties of clusters and of the galaxies they 
contain. 

Briefly, WINGS is a long term multiwavelength project 
based on deep optical (B,V) wide field images (~ 35'x35') of 
77 fields centered on nearby clusters of galaxies selected from 
three X-ray flux limit ed samples compiled from R OSAT All-Sky 
Surve y data ( BCS. lEbeling et alj 119961: eBCS, lEbeling et all 
119981: XBACs. libeling et al.ll2000b . The selected redshift range 
(0.04 < z < 0.07) was set to balance a wide linear field of 
view (~ 1 .4 Mpc x 1 .4 Mpc at z ~ 0.04) and a high spatial 
resolution (1.34kpc/" at z ~ 0.07). To reduce the effects of 
galactic extinction in our analysis, only clusters located far 
from the galactic equator (\b\ > 20°) were kept. The obser- 
vations were carried out using the wide field cameras of two 
telescopes: the 2.5m Isaac Newton Telescope (WFC@INT) 
and the MPG/ESO-2.2m telescope (WFI@ESO). The op- 
tical data were complemented wi th spectroscopic follow 
up of a subsample of 48 clusters dCava et all 120081) using 
WYFFOS@WHT (A range=3800-7000 A, A resolution=3A) 
and 2dF@AAT {A range=3600-8000 A, A resolution=6A). 
Also, we obtained J and K imaging of 32 WING S clusters 
with WFC AM @ UKIRT dValentinuzzi et all 12008b and we 
are collecting U/H band imaging of WINGS clusters with 
wide-field cameras at different telescopes (INT, LBT, Bok), that 
will be used to analyse the stellar masses and of star formation. 

iPaper J started the series devoted to the analysis of the optical 
data within the so called WINGS-OPT subproject. It presented 
the criteria followed to construct the initial sample of clusters 
and described the optical observations (B and V bands) and their 
photometric and astrometric quality. 

This is the second paper in the series regarding the optical 
photometry of the WINGS project and it is devoted to the re- 
lease of the catalogs containing the basic photometric parame- 
ters of all the objects found in the fields of 77 nearby clusters of 
galaxies. Also, it describes the procedures that were followed to 
construct the catalogs (Section [2). The process includes a pre- 
liminary treatment of the images (§ 12. U that has allowed us to 
improve the photometry of many galaxies, especially the largest 
ones. In Section[3] catalogs are presented and their overall qual- 
ity is checked including their photometric quality (§ 13.2. U . com- 
pleteness and star/galaxy separation (§ 13.2.31 ). A final summary 
will be found in Section[4] 

Throughout this paper we will use a cosmological model 
with parameters: Ho = 75 kms~' Mpc -1 , Q.m = 0.3 and Q.\ = 
0.70. 



1 A more det ailed comparison from the spectroscopic point of view 
will be made in lCava et al]<2008ft . 

2 Be aware that in IPaper 5 we used instead Hq = 70 kins' Mpc -1 , 
however this change doesn't modify the results showed there 
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Fig. 1. Example of the problems of SExtractor to produce the correct segmentation of the extended galaxies. From left to right: (a) 
BCG of cluster A193 in the V band image; (b) segmentation performed by SExtractor without subtracting the galaxy; each color 
represents those pixels assigned by SExtractor to different objects ; (c) segmentation of the remaining objects after subtraction of 
the BCG; (d) final segmentation of the BCG. 



2. Detection, basic photometry and star-galaxy 
classification 

The whole process of source detection, computation of their 
basic photometric parameters and star-galaxy classificatio n 
was performed using SExtractor dBertin & Arnoutsl Il996l) . 
SExtractor is a very efficient tool to find sources in large for- 
mat images. Among its characteristics the most interesting are its 
ability to separate blended sources and to estimate the probabil- 
ity of an object of being a star or an extended source. However, 
both processes need to be fine-tuned by input parameters that 
are fixed for each single image. To improve the final outcome 
of SExtractor a preliminary treatment of the images was done 
before running the source detection program as well as a final in- 
teractive checking of the star-galaxy classification, as described 
in the following. 

2.1. Preliminary image treatment 

In images like ours in which the ranges of size and brightness 
are quite large, it is impossible in practice to find a single set of 
values of SExtractor's input parameters that could work ade- 
quately for all the objects. The extreme situation is when a large 
galaxy is contaminated by a number of small projected sources 
as is the case of the brightest central galaxies (BCGs) of our sam- 
ple of clusters. This issue is illustrated in FigQ]with the BCG 
of the cluster A193 (panel (a)). The segmentation^] of the cen- 
tral galaxy and the objects projected onto it is shown in panel 
(b). Clearly, the program has erroneously assigned pixels of the 
large galaxy to the small projected objects producing poor pho- 
tometry of the large galaxy as well as of the objects projected 
onto it. Something similar happens with very bright stars whose 
extended halos affect the photometry of the nearby objects. 

To minimize the effect of such large halos in the photometry 
of close objects as well as to improve the photometry of the ex- 
tended galaxies themselves we have developed a custom-made 
procedure. It consists of first the modeling and removal of these 
halos before running SExtractor on the image. This improves 
the photometry of the small projected objects. In a second step, 
an image is constructed containing only the extended galaxies 
and the photometry is performed on it with SExtractor. In this 
last image the pixels of the projected small objects are replaced 
by the values of the models and the rest of the pixels are left 
unchanged. This replacement reduces the contamination of the 
projected objects in the photometry of the large galaxies. An ad- 
ditional advantage of the procedure is that the removal of the 

3 The segmentation is the way in which a program of source detection 
assigns pixels to each object. 



extended halos also improves the determination of the global 
background map. 

In the following we give a brief description of the process 
that is explained in more detail in AppendixlAl 

The process starts with the computation of a first background 
map which is subtracted from the image. Then, the problematic 
objects (galaxies and stars with extended halos) are located. For 
each of these objects, a mask of theprojected objects is con- 
structed and elliptical isophotes are fiQ The resulting fit is used 
to construct a model of the halo that is afterwards subtracted 
from the original image (i.e. the image before the background 
subtraction). After doing that with all the selected objects, the 
resulting image without the large halos is used to refine the back- 
ground map as well as the masks of the smaller objects. With the 
new background map and masks the process is repeated to im- 
prove the final results. It is found that one iteration is enough to 
achieve sufficiently accurate photometry. 

The last step is the construction of a complementary image 
containing only the previously removed galaxies (of course, the 
bright stars are also avoided in this image). It is important to 
note that the photometry of the large galaxies is not performed 
on the models but rather using the original pixels. The models 
are used only those regions occupied by projected galaxies or by 
interchip regions. 

At the end of the whole process, which is run in V band as 
well as in B band, the output is two background subtracted im- 
ages in each band. One of the images contains all the objects ex- 
cept the largest ones and the other one only the extended galaxies 
that were removed from the first image. 

The improvement of the procedure is illustrated in the last 
two panels of FigQ]in which we show the segmentations of the 
projected objects (panel (c)) and of the BCG (panel (d)) after 
applying our procedure. 

2.2. Detection and basic photometry 

To construct the photometric catalogs, the detection of sources 
and the basic photometry relies on SExtractor and more 
specifically on the modified version SExtractor 2.3.2 by G. 
MorrisorQ following the changes initially made by B. Holwerda 
for SExtractor v.2.2.2.These versions have the advantage of 
computing additional parameters related to the light distribution 



4 To achieve a successful fit a preliminary 'masking' procedure 
is needed of the objects that appear projected onto the halos. See 
AppendixlAlfor more details. 

5 http://www.cfht.hawaii.edu/~morrison/home/SExtractor.html 
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Fig. 2. Four different kinds of diagrams which can be used to 
\og(Area) ; log(FWHM) - log(Area); V - (V(r < 1") - V). The 
A2589 (FWHM, = 0."85; upper panels) and A2626 (FWHM* = 1." 
green dots unclassified objects. For better visualization, only 20 
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perform a star/galaxy separation. From left to right: fi max -V ; V- 
data come from two fields with quite different seeing conditions: 
64; lower panels). Red squares are galaxies, blue crosses stars and 
% of the objects have been plotted. 



in the obje cts such as the concentratio n, the contrast and the 
asymmetry dAbraham et al.Lll994l[T996b . 

The detection of sources, the determination of their posi- 
tions and geometrical parameters as well as the star/galaxy clas- 
sification was done on the V band image. Since we had im- 
ages from two different cameras with different pixel scales, 
instead of setting a single detection threshold per pixel for 
both cameras we opted for setting the same detection thresh- 
old per square arcsec. The advantage of this approach is that 
with similar observational conditions the limiting surface bright- 
ness is independent of the pixel size. The disadvantage is that 
in the images with lower detection threshold per pixel (in our 
case the WFI@ESO) the rate of spurious detections at low 
brightness is higher. However, since we were interested in the 
bright part of the luminosity function, where the galaxies of 
the cluster outnumber those of the background, it was desir- 
able to perform the photometry of the galaxies within similar 
surface brightness limits. In order to balance photometric depth 
(i.e. detection of low surface brightness galaxies) with a low 
rate of spurious detections, the detection threshold was set to 
4.5 cr^/arcsec 2 (o> g is the standard deviation of the background 
signal) which corresponds to 1.5 cr^/pixel for the WFC@INT 
and 1.07 cr/,g/pixel for the WFI@ESO. Given the typical values 
of (Tb g found in our images these limits translate in a detection 
limit of HThreshoidY) ~ 25.7 mag/arcsec 2 . As a comparison, at 
the same signal to noise level the images from the SDSS0 in 
g band reach finreshoidig) < 25.2 mag/arcsec 2 while those in r 
band reach iijhreshou(f) < 24.7 mag/arcsec 2 . 

SExtractor was run twice with the B band images using 
both the single-image and the dual-image modefl The cata- 
log resulting from running SExtractor in single-image mode 



6 The following values have been calculated from the image identi- 
fied by (run,rerun,camcol,field)=( 1889,40,3, 105) which is the one taken 
in most favorable conditions 

7 In single-image mode the detection and the photometry are done 
using the same image. In dual-image mode the detection is done in one 



allowed us to reduce the number of spurious detection be- 
cause only objects with detections in both bands (separated 
by less than 1767 = 5 pixel (WFC@INT) = 7 pixel (WFI@ESO)) 
were kept. However, to obtain correct color indexes of the galax- 
ies we performed a second run of SExtractor on the B band 
image using the dual-image mode in which the V band image 
was set as the detection image. This procedure ensured that the 
measurements of each galaxy in both filters were done in the 
same regions. 

SExtractor was run without using weighted images be- 
cause preliminary tests showed that the photometry was more 
robust when not using weighting images^ As a consequence, 
photometric errors computed by SExtractor are unrealisti- 
cally small. The estimation of these errors is explained in 
Section 

To avoid the influence of the interchip regions, the pixels in 
these regions were set to zero and the objects falling within a 
band of several pixels to one of the edges of the chips^ were 
excluded from our catalogs. 

The saturated star£3 were not included in the catalogs be- 
cause of the uncertainty in their photometry. 

Once having obtained the photometry in both bands, we 
applied the color term correction of the flux calibration. This 
was done using the color index (B-V) measured in apertures of 



reference image and the photometry of the detected objects is done in a 
second image. 

8 Nowadays, the use of weighted images is widely used and a safe 
procedure, however, at the moment of constructing the catalogs we saw 
that the use of the weighted images increased the uncertainties in the 
magnitudes so we opted for not using them. 

9 The bands were constructed interactively to ensure that the prob- 
lematic regions were avoided. 

10 The exposure time of the single exposures were set to avoid the 
saturation level for galaxies so no galaxy was found to reach such a 
level. 
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R=5 kpc, except for objects with smaller linear radius for which 
smaller apertures (R=2 kpc) were used. 

For each field, SExtractor was run on the images obtained 
for both the bright and faint galaxies (see Sect |2. lb and the final 
catalog was obtained by merging the two catalogs. 

2.3. Star-galaxy classification 

The star-galaxy classification was made relying upon the stel- 
larity index (CLASS_STAR) computed by SExtractor which 
ranges from (galaxies) to 1 (stars). We decided to make a 
rather robust selection in two groups, "stars" and "galaxies", with 
the highest probability of being correct. Objects not fitting those 
criteria were left in a third group of "unknown" classification. 
Practically we imposed the following criteria: 

Stars CLASS.STAR > 0.8 

Galaxies CLASS.STAR < 0.2 
Unknown 0.2 < CLASS_STAR < 0.8 

In spite of the robustness of the criteria that we used, the fi- 
nal catalogs were checked for possible star/galaxy misclassifica- 
tions. Since the number of objects was too large for any individ- 
ual analysis, we opted for using several plots of different combi- 
nations of parameters to detect these possible misclassifications. 
In Fig. [2] we show some of the combinations of parameters that 
separate stars from galaxies: 

- V-p malc , where fi max is the surface brightness of the brightest 
pixel in an object and V is the total V band apparent magni- 
tude which actually is SExtractor's MAGJVUTO. 

- log I0 (7 sophotal Area) - V, where Isophotal Area is the area 
in pixels of each object above the threshold. 

- log 10 (/ sophotal Area) - log(FWHM). 

- V - (V(r < 1") - V), where V(r < 1") is the magnitude 
measured in an aperture of radius=l". 

In all these diagrams, stars populate a narrow and well de- 
fined region while galaxies are more spread throughout the 
plane. This clear segregation decreases for smaller and fainter 
objects for which stars and galaxies populate similar regions in 
the diagrams. Therefore, these diagrams are not really useful for 
faint (V > 20) and small objects (log{Area(pixels)) < 2). 

Displaying these kinds of diagrams for each cluster field it 
was easy to detect outliers that most probably represented ob- 
jects with a wrong classification as well as remaining spurious 
detections. The (very few) suspicious objects were then visually 
checked and their classification changed if needed. 

At the end of this process we computed the degree of mis- 
classification of the original SExtractorcatalogs. It was found 
that the fraction of misclassification was minimal, less than 1% 
of misclassified stars and < 0.6% of misclassified galaxies, up 
to V ~ 22. Therefore, we consider that the remaining misclassi- 
fications in this range of magnitude should be even smaller. For 
fainter objects the previous diagrams lose their utility since stars 
and galaxies populate almost the same region and even the visual 
check becomes useless. 

The published catalogs will be regularly updated to 
correct for possible new spurious objects or misclassifi- 
cations found so users are encouraged to use the lat- 
est version that will be published in WINGS web site: 
http : //web . oapd. inaf . it/wings/. 



3. The catalogs 

3.1. Catalogs Description 

We have constructed three catalogs for each field containing re- 
spectively galaxies, no saturated stars and objects of unknown 
classification. The structure of all of them is the same and an 
example is shown in Table Q] 

The parameters stored for each object are the following (in 
parentheses we give the name of the output parameter from 
SExtractor used to calculate them): 

- Equatorial coordinates (J2000.0) of the barycenter 
(X_IMAGE , Y.IMAGE) and of the peak of emission 
(XPEAKJMAGE , YPEAKJEMAGE), taken from the V band 
image. 

- Isophotal area (ISOAREAJMAGE). 

- Kron radius (KRON_RADIUS). 

- Full width at half maximum (FWHPLIMAGE). 

- Axis ratio computed as the ratio of the SExtractor's pa- 
rameters B_IMAGE and A.IMAGE. 

- Position angle with respect to the North and measured 
counter-clockwise (THETA.IMAGE). 

- SExtractor's FLAG 

- SExtractor's stellarity index (CLASS_STAR) 

- Concentration index measured as the ratio between the 
flux in the central 9%0 of the pixels and the total flux 
(CONCENTRATION). 

- Contrast index measured as the ratio between the flux in the 
brightest 30% of the pixels and the total flux (CONTRAST). 

- Surface brightness of the brightest pixel (MUJ1AX). 

- And the following magnitudes in both bands: 

- Total magnitudes: SExtractor's MAG.ISO, MAG_IS0C0R 
and MAG_AUT0 

- Three magnitudes at fixed physical apertures (at the tar- 
get cluster's redshift): R=2kpc, 5 kpc, 10 kpc. 

- Three magnitudes at fixed angular apertures: l'/6, 2'.'0, 

The photometric catalogs are public only in elec- 
tronic format via CDS and at WINGS' website 
http://web.oapd.inaf.it/wings/. The reduced im- 
ages in both bands will be made available upon request to the 
authors. 

In the following we give indications about some points to be 
taken into account when using the catalogs. 

Coordinates The a strometry of the images is discussed in 
Appendix A.4 of iPaper ]Pl . The overall quality is very sat- 
isfactory (see Appendix A.4 of Paper I) with uncertain- 
ties of the order of 0703 for stars and ranging from 0718 
for WFI@ESO images to 0725 for WFC@INT ones. For 
each object two positions have been computed: that of the 
barycenter of the light distribution and that of the peak of 
the emission. For most objects, they are almost coincident. 
However there are two situations in which it is better to use 
one and not the other. First, for faint low surface brightness 

11 Due to a bug found in the modified version of SExtractor used to 
compute the concentration index, the flux is not measured in the central 
30% of the pixels but in the central 9%. 

12 These fixed apertures have been chosen for the multifiber spec- 
troscopy of the WINGS fields. 176 is the projected diameter of the fibers 
in Autofib2@WHT while in the 2dF@AAT the diameter varies radially 
in the field from 27 16 in the center to 270 in the edges. 

13 The astrometry was done with respect to the USNO-A2 Catalogs 
which result in slight differences (< 0"3) with respect the USNO-B1. 
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Table 1. Example of entries in one of the photometric catalogs. 



ID 


<*Bary 




"Peak 


Speak 


Area 


TKron 


FWHM 


b/a 


PA 


FLAG 


S.I. 




(deg) 


(deg) 


(deg) 


(deg) 


(arcsec 2 ) 


(") 


(") 




(deg) 






WINGS J004135.4-090100.5 


10.39751 


-9.01681 


10.39751 


-9.01682 


17.39 


1.49 


6.26 


0.73 


63 


19 


0.02 


WINGS 1004106.0-090104.2 


10.27495 


-9.01783 


10.27500 


-9.01781 


41.86 


0.83 


2.03 


0.75 


31 





0.03 



Vband 


Cone. 


Contr. 




MAG.IS0 


MAG.ISOCOR 


MAGJUJTO 


MAG(2Kpc) 


MAG(5Kpc) 


MAG(lSKpc) 


MAG (17 6) 


MAG (278) 


MAG (27 16) 


0.384 


0.543 


22.97 


20.71 


20.44 


19.86 


21.43 


20.63 


20.10 


22.53 


22.19 


22.09 


0.460 


0.777 


21.03 


18.58 


18.55 


18.57 


19.11 


18.60 


18.51 


20.48 


20.09 


19.97 



B band 


MAG.IS0 


MAG.IS0C0R 


MAG_AUT0 


MAG(2Kpc) 


MAG(5Kpc) 


MAG(18Kpc) 


MAG (176) 


MAG(278) 


MAG (27 16) 


21.41 


20.46 


20.00 


22.18 


21.09 


20.21 


23.37 


23.01 


22.90 


19.60 


19.48 


19.58 


20.16 


19.61 


19.47 


21.58 


21.18 


21.05 



ID: Object internal identification. 

{(^Bary^Bary) '■ Equatorial coordinates (J2000.0) of the barycenter. 
(piBeay^Bary) '■ Equatorial coordinates (J2000.0) of the brightest pixel. 
Area: Area above the detection threshold. 
r Kmn '■ Kron radius used to computed the MAG_AUT0 magnitude. 
FWHM: Full width at half maximum assuming a Gaussian core, 
b/a : Axis ratio. 

PA : Position angle of the major axis (North=0°, counter-clockwise). 

FLAG : SExtractor's FLAG 

Si. : SExtractor's stellarity index CLASS.STAR 

/i mflr : Surface brightness of the brightest pixel. 

Cone. : Abraham's concentration index, measured as the ratio between the flux in the central 9% of the pixels and the total flux. 
Contr. : Abraham's contrast index, measured as the ratio between the flux in the brightest 30% of the pixels and the total flux. 
MAG_IS0 : SExtractor's isophotal magnitude. 
MAGJSOCOR : SExtractor's isophotal corrected magnitude. 
MAG_AUT0 : SExtractor's Kron aperture magnitude. 

MAG(2Kpc) , MAG(5Kpc) , MAG(lQKpc) : Magnitudes within apertures of radius 2Kpc, 5Kpc and lOKpc, respectively, measured at the clusters' redshift. 
MAG (176) , MAG (278) , MAG (27 16) : Magnitudes measured within apertures of diameter equal to 176, 270, 27 16, respectively. 



galaxies the barycenter position is preferable since the posi- 
tion of the peak is highly affected by the noise. On the other 
hand and most importantly, for large extended galaxies (such 
as the brightest cluster galaxies) the coordinates of the peak 
emission are a better choice to locate the center of the object 
because the barycenter is strongly affected by the shape of 
the most external isophotes. 

Shape parameters Care must be taken when using the shape 
parameters especially for objects in crowded fields and/or 
contaminated by close companions. The use of the FLAGS 
values can help to check this. 

Total magnitudes SExtractor gives four different measures 
of total magnitude. MAG_IS0 is the integrated light above 
the detection threshold. MAG_IS0C0R is an attempt to cor- 
rect for the light lost in the wings of the objects mak- 
ing use of a Gaussian model approximation. It was found 
that for B band magnitudes computed in dual mode 
MAG_IS0C0R can give wrong values so we discourage its us- 
age. MAGJUJTO is another attempt to guess t he tota l mag- 
nitude of the galaxies using Kron aper tures (iRronl Il980t 
llnfantel Il987t iBertin & Arnoutsl 119961) . However, several 
works (e.g. iFranceschini et all Il998h have shown the prob- 
lems of this approximation, especially its dependence on the 
light profile of the object. In Section [3.2.1| we will show that 
while for stellar and exponential profiles the MAGJUJTO mag- 
nitude is a good approximation, for de Vaucouleurs profiles 
there is a systematic offset of ~ O.2"0 Finally, SExtractor 
includes another total magnitude called MAGJ3EST magni- 
tude. This corresponds to MAGJUJTO except in those cases 



in which the light of the object has a contamination by ex- 
ternal sources higher than 10%, when MAG.ISOCOR is used 
instead. Several sources including the last SExtractor's 
Users Guide discourage the usage of MAG_BEST in favor of 
MAG_AUT0. Our experience when constructing the catalogs 
has convinced us that MAG_AUT0 gives the best approxima- 
tion to the total magnitude of objects. 
Saturated stars Since the aim of the catalogs is the study of 
clusters of galaxies, the saturated stars were discarded. This 
should be taken into account when using our catalog of stars. 
In the case of WFC@INT fields the incompleteness pro- 
duced by saturation affects stars brighter than V ~ 17, while 
for WFI@ESO fields this limits is V ~ 16, the difference 
mainly from the different pixel size. 



In addition to the information on the single objects found in 
each field, in AppendixlDlthe reader can find a set of tables with 
additional information about the target clusters, peculiarities of 
the fields and the conditions of observation that could be useful 
when working with our published data. 



3.2. Quality checks 

In this section we describe the internal photometric qualitjH of 
the data and the accuracy of the star/galaxy classification in our 
catalogs. 



14 A detai led explanation of t he Kr on magnitude and its problems can 
be found in lGraham & Drivej d2005h 



15 The comparison with external data (SDSS) was already done in 
iPaper J. section 5.4. 
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Fig. 3. Photometric errors computed from simulations. Upper 
panel: Systematic offset of the measured magnitude (V out p U t) 
with respect to the input magnitude {Vi„ put ) as a function of 
the input magnitude. Lower panel: Dispersion of this difference. 
Green squares show the differences between the photometry 
of two fields (A780 and A970) observed with W FC@INT and 
WFI@ ESO (same as lower fifth panel of Fig.l 1 in lFasano et alJ 
(120061) . 



3.2.1. Photometric errors 

When running without weighted images, SExtractor computes 
just the photon-noise errors of the magnitudes. For our deep im- 
ages these values are unrealistically small and for this reason 
we do not report such errors in our catalogs. In fact, the er- 
rors in the photometry are mostly related to local background 
variations or to contamination by the light of close compan- 
ions. To estimate these uncertainties we performed simulations 
in which synthetic stars as well as galaxies with exponential 
and de Vaucouleurs profiles were inserted in the original im- 
ages. Fig. [3] (upper panel) shows the differences found between 
the input magnitude (Vi„ put ) and that measured by SExtractor 
(V ma put) as a function of the input magnitude for the three 
types of objects used in the simulations: stars and the two types 
of galaxy models. Each point is th e central biweight estima- 
tor lliirelFiynnXGibhardl [l990|) of the difference {V inpul - 
y output) for objects in an interval of A(Vi„ p „i) = 0.3'". In the lower 



panel we show the scale biweight estimator of V, 



input 



V, 



Hit pill 



m 



the same intervals. Note the use of a logarithmic scale in the 
ordinate axis. 

It is clear that for de Vaucouleurs profiles, SExtractor pro- 
duces total magnitudes that are ~ 0.2'" fainter than real ones, or 
equivalently, fluxes tha t are ~ 20% lower. This is a known prob- 
lem o f SExtractor dFranceschini et all Il998t iBenftez et al.L 
I20041) although the differences reported between the real and 
the measured magnitudes vary from work to work. Contrary to 
previous works, we find that the difference between input and 
output magnitudes remains almost constant up to our detection 



limit (V ~ 23). On the other hand, stellar and exponential pro- 
files show good agreement with mean differences < 0.05'". 

The simulations also give us a more realistic measure of the 
uncertainties on the magnitudes (lower panel of Fig . [3} . At fixed 
input magnitude, the largest uncertainties are for galaxies with 
de Vaucouleurs profiles and the lowest are for stars. The trend 
with input magnitude is similar for the three types of objects with 
a steeper increase of the uncertainties at magnitudes fainter than 
V ~ 20. As a comparison the figure also shows (green points) the 
differences in magnitudes found when comparing the data from 
two fields (A780 and A97 0) observed w i th our two instrumental 
set-ups (see Fig. 11 from iFasano et a il (120061) ). Therefore, the 
typical uncertainties in our catalogs up to V ~ 21.5 is less than 
0.1 mag. For magnitudes fainter than V = 21.5, it is possible to 
estimate the typical uncertainties with the following expressions: 



cr v (S tar) 
cr v (Exp) 
<Tv(deVauc) 



10 
10 
10 



,0.309V-7.62 
,0.271^-6.67 



0.259V-6.34 
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Fig. 4. Comparison of the values of the total magnitude (lower 
panels), color index (B - V)(middle panels) and isophotal area 
(upper panels) measured by SExtractor of the modeled galax- 
ies before and after separating them from the rest of the ob- 
jects. Left: Distribution of the differences of these parameters 
as a function of the total V magnitude. Right: Distribution of the 
same differences. Index "B" denotes the values before applying 
our procedure and index "A" indicates the values after applying 
our procedure. 



3.2.2. Photometry of the modeled objects 

In this section we check the improvement in the photometry of 
the objects modeled after applying the procedure described in 
AppendixlAlto the images. 

Figure|4]shows the effect of the model subtraction on the out- 
put from SExtractor of some photometric parameters of the 
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Fig. 5. Average detection rate in each observing run computed 
from simulations. The stronger black line is the detection rate 
averaged over the 77 fields. 



modeled galaxies. The large increase in the isophotal area after 
applying our procedure clearly reflects the issue with the seg- 
mentation of large galaxies (see Fig. QJ>). This also affects the 
computation of the total magnitudes, as can be seen in the lower 
panel of Fig. [4] On the other hand, the color index (B-V) shows 
smaller variations because of the use of aperture magnitudes to 
compute it. 

In addition to the improvement in the photometry of the sub- 
tracted objects we also increased the detection rate of objects 
projected onto the modeled halos. In the case of objects pro- 
jected onto the halos of the BCGs we found ~ 16% more objects 
after the subtraction of the BCG. This is a relevant result es- 
pecially in studies of the spatial distribution of galaxies in the 
cluster. 

3.2.3. Completeness and goodness of the star/galaxy 
classification 

The simulations described above were originally intended also to 
check the detection rate (completeness) and the success rate of 
the star/galaxy classification. However, we realized that the sec- 
ond step could not be done due to an unexpected issue. We found 
that commonly used programs to simulated objects (stars and 
galaxies) do not work well at faint magnitudes. The result is that 
when running SExtractor, simulated stars have a lower prob- 
ability of being well classified than real stars (see Appendix IBl. 
For that reason, simulations were used only to estimate the de- 
tection rate in each field. 

Figure[5]shows the average detection rate run by run (dotted 
curves) and the total average detection rate (continuous line). 
The curves are quite similar and we can considered our global 
catalogs 90% complete for objects with V < 21.7. The 50% 
detection level is reached at V ~ 23.2. In addition, we provide 
in Table |D.2l the V band magnitudes at which the detection rate 
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Fig. 6. Spatial variation of the stellar FWHM in the differ- 
ent runs of WINGS. Each point is the cent ral biweight esti- 
mator (CBE jBeers. Flvnn & Gebhardl d 19901) ) of the ratio be- 
tween the FWHM (FWHM*) and the FWHM used as input 
in SExtractor (FWHM SEk ) for stars in bins of 500 pix- 
els (WFC@INT) or 700 pixels (WFI@ESO) along both axis. 
The errors have been computed as the scale biweight estimator 
(SBE). 



drops to 90%, 75% and 50% in each field as well as the surface 
brightness detection thresholds in the same band. 

As an external check of the completeness computed with the 
simulations we took the image of the center of A 1795 taken from 
the Space Telescope and compared the detection rates in that 
image with ours. The results are presented in Appendix ICl and 
they are in good agreement with those from the simulations. 

Concerning the star/galaxy classification, for objects with 
V < 21 the difference between light profiles of stars and galax- 
ies is large enough to allow SExtractor to distinguish them 
easily and, moreover, the diagrams like those in Fig. [2] can be 
used to reclassify the cases where SExtractor fails. Therefore, 
for these bright objects the number of misclassifications can be 
considered negligible. 

The problem arises for fainter objects whose profiles be- 
come more similar and small variations in the local conditions 
(e.g. small variations in the background) can make SExtractor 
fail to perform a good classification. Moreover, the number of 
objects is so large that interactive cleaning is unpractical and 
even the visual check cannot help. This is worsened by the fact 
that the point spread function actually varies throughout the field 
(Fig. [6]). For these reasons, a statistical approach appears to be 
the best option to check the reliability of the star/galaxy classifi- 
cation at faint magnitudes. As stated above, the simulations were 
found to be unreliable at faint magnitudes so we opted for an- 
other way to estimate the goodness of our star/galaxy classifica- 
tion as well as to guess the composition of the "unknown" class 
of objects. First of all, we can compare the fraction of objects in 
each group (stars/galaxies/unknown) as a function of the V mag- 
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Fig. 7. Relative numbers of the three types of objects: galaxies, 
stars and unknown objects. WFI@ESO fields show a higher frac- 
tion of stars because some of them point closer to the direction 
of the galactic center. 



nitude (see FigO. As expected, the fraction of unknown objects 
increases at fainter magnitudes. This increase rises more steeply 
for V > 22, i.e. when the completeness starts to fall rapidly. A 
interesting point is that this rise corresponds to a change in the 
trend of the fraction of galaxies, pointing to the fact that most 
unknown objects should be galaxies. Another way to see this 
is by relying upon external sources for galaxy and star counts. 
For the galaxies we took d ata from the ESO- Spitzer Imaging ex- 
tragalactic Survey (ESIS, iBerta et ail 120061) while for the stars 
we used the models of distri bution of stars in th e Galaxy of the 
Observatory of Besancor P'l dRobin et all [2003b . Of course, this 
check is relevant only in the faintest part of the magnitude dis- 
tribution where counts from stars and background galaxies dom- 
inate those of the target cluster. With both counts we computed 
the fraction of galaxies with respect to the total number of ob- 
jects as a function of the apparent magnitude (continuous line in 
Fig- Si. This is compared with the same fraction measured from 
WINGS fields (dashed lines) although measured in off-center re- 
gions to avoid the high contribution of galaxies from the cen- 
ter of the clusters. If the "unknown" objects are ignored (short- 
dashed line in the figure) the fraction of galaxies observed in our 
fields becomes smaller than that observed in the ESIS+Besancon 
sample for V > 20.5. This is an unexpected result since our fields 
include the contribution of the cluster. To check the influence of 
the classification of galaxies as "unknown" objects, we assumed 
that all the "unknown" objects were galaxies and recomputed the 
fraction (long-dash line). This produces an overestimation of the 
fraction of faint galaxies. It is possible to compute the minimum 
fraction of galaxies in the "unknown" sample which makes the 
WINGS galaxy fraction equal that of the background counts of 




http://bison.obs-besancon.fr/modele/ 



Fig. 8. Lower panel: Average fraction of galaxies as a function 
of the V magnitude. Continuous line: Fraction computed assum- 
ing the stell ar counts from the models of Besancon and galaxy 
counts from lBerta et al.ld2006l) . Short dashed line: Data from off- 
center regions ignoring the unknown objects. Long dashed line: 
Data from off-center regions including the unknown objects as 
galaxies. Upper panel: Minimum fraction of galaxies in the un- 
known objects. 



galaxies and stars. This is shown in the upper panel of Fig.[8]and 
for V > 23 this fraction is greater than 80%. 

The agreement between the fraction of galaxies observed and 
expected from external sources reassures us that even at very 
faint magnitudes (V ~ 24) the star/galaxy classification is not a 
mere random assignment from which one would expected half 
stars and half galaxies. Since the difference in the fraction of 
galaxies is always < 10% between our data and those from the 
external sources we can assume that statistically the star/galaxy 
classification holds quite well up to V ~ 24. 



4. Summary 

In this paper we present the first data release from the WINGS 
project. We have produced photometric catalogs of all the ob- 
jects found in 77 fields centered on nearby clusters. The objects 
have been classified as galaxies, stars and objects of unknown 
type. For each object we give positions, geometrical parameters 
and different kinds of total and aperture magnitudes in both ob- 
served bands, B and V. 

Our catalogs are 90% complete at V ~ 21.7 and 50% at 
V ~ 23.2 although these values vary from field to field, espe- 
cially depending on the total exposure time. Our star/galaxy clas- 
sification relies on the CLASS.STAR parameter of SExtractor 
but we have performed a visual checking of the brightest ob- 
jects (V < 22) that makes the misclassifications negligible. For 
the faintest objects we tried to use simulations but we found that 
these were not reliable in the appropiate range of brightness (V > 
22). However, when we compare the fraction of galaxies found 
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in our catalogs with those from external sources (ESIS' deep 
galaxy counts and stars counts from Observatory of Besancon's 
Galaxy Model), we find that both are compatible within +10%. 
This result makes us confident in our star/galaxy classification 
even at faint magnitudes, at least statistically. Again, the relia- 
bility of the classification varies from field to field, depending 
on the seeing and also on the spatial variation of the PSF. 

A notable feature of these catalogs is that extended galax- 
ies have been treated in a special manner which has allowed us 
to improve their photometry, with several cases in which the to- 
tal magnitudes were initially underestimated by more than one 
magnitude. This is a critical issue for our fields, since our clus- 
ters are often centered on very bright, extended galaxies, which 
are themselves important targets of our study. We have been able 
not only to improve the photometry of these galaxies but also to 
detect ~ 16% more objects around the BCG relative to the case 
in which the galaxy is not subtracted. 

In addition to the catalogs we also release the reduced im- 
ages in both bands (B,V) and the preprocessed images with the 
large halos removed as well as the images with only the largest 
galaxies. 

The scientific analysis of these data will be published in the 
subsequent papers of this series. 
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Appendix A: Full description of the extended halo 
removal 

In this appendix we give a detailed description of the whole pro- 
cedure followed to remove the objects with extended halos. 

This is the outline of the whole process, while the single 
steps are explained in some detail thereafter: 

1. A background map is computed using the original images. 

2. The extended objects (galaxies and stars) to be removed are 
selected interactively. 

3. For each selected object, the following steps are performed: 

(a) The objects projected onto it are maskecQ with an inter- 
active procedure. 

(b) Its isophotes are fitted with the IRAF task ellipse. 

(c) The resulting elliptical isophotes are then used as the in- 
put of the task bmodel to construct a model of the ex- 
tended object. 

(d) The model is subtracted from the original image. 

(e) In the case of galaxies, when the subtraction leaves resid- 
uals due to the structure of the galaxy, these are removed 
manually with imedit. 

4. After removing all the selected objects from the initial image 
a new background map is calculated. 

5. A mask of the remaining objects is constructed with 
SExtractor. This mask is then used in a second iteration 
to improve the interactively-made masks (point fl3a| l) applied 
during the isophote fitting procedure. 

6. At this point, the procedure can be repeated from point (0 
until the subtraction is satisfactory. 

The process is iterative and could be repeated as many times 
as necessary but we found that one iteration was enough to reach 
a good photometric quality. 

Now we explain in more detail all the steps of the procedure: 

First background computation Due to the variations found 
from chip to chip and to avoid the influence of the interchip 



17 The masking procedure consists of creating an auxiliary image or 
mask that indicates to the fitting program which pixels avoid to be in the 
fitting procedure. This can be done, for example, by setting to zero all 
the pixels of the mask except those whose position correspond to pixels 
in the original image that should not be used in the fit. 
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Fig. A.l. Upper panel: Fraction of galaxies that were modeled. 
Lower panel: Total number of galaxies that were modeled. 



regions (which SExtractor treats as if they were part of the 
actual image) the images were split into their original chips. 
After that, an initial background estimation was computed 
for each chip using SExtractor with BACK_SIZE=256 and 
BACK_FILTERSIZE=3 as input parameter values. To reduce 
the influence of bright pixels, those with intensities above an 
established threshold were replaced by the mode of the pixel 
intensity of the same chip, which can be considered as a first 
estimation of the background level. The threshold was set 
manually before starting the process seeking to not remove 
the bright part of the background. 

Once the background was obtained for each chip, the sin- 
gle background images were mosaiced to construct a global 
background map. This background map was then subtracted 
from the original image. 
Selection of the objects to be removed The selection of the 
objects to be modeled and removed was done in a subjec- 
tive way but following several guidelines. First, the BCG was 
always modeled. Elliptical galaxie{3 highly blended and/or 
surrounded by small objects were also modeled. The selec- 
tion of the stars to be modeled was more field-dependent. 
As the amount of work needed increased substantially with 
the number of objects to be modeled, in those fields with a 
high density of stars only the brightest stars were modeled. 
Sometimes, a star was modeled to avoid the contamination 
from a close bright galaxy. 

Fig. lA.ll shows the distribution in V of the galaxies that have 
been processed in this way. Up to V ~ 14 most of the galax- 
ies have been modeled. 
Fitting and modeling To reduce the computation time, the pro- 
cess of isophote fitting, modeling and model subtraction was 



18 Since the process is based on the construction of models of elliptical 
isophotes it is not suitable for spiral or irregular galaxies. 



done on small images of the selected objects extracted from 
the background-subtracted global image. 
The isophote fitting was done using the IRAF task ellipse. 
However, before that, all the projected objects and problem- 
atic regions (e.g. interchip regions) were masked. Due to the 
large number of objects that should be masked and taking ad- 
vantage of the iterative process, this first step was enough to 
make a rough mask. Below we will show how the improve- 
ment is achieved. 

A few differences were introduced when dealing with stars 
or galaxies. For stars the ellipticity was fixed to zero and the 
fitting did not reach the innermost region (usually saturated 
and, therefore, not suitable for fitting isophotes). Also, the 
center of the isophotes was not fixed because quite often the 
reflections of the stars were off-center with respect to the 
central regions. For galaxies no restrictions were imposed. 
The output of ellipse was used as the input of the task 
bmodel which allows us to construct a two dimensional 
model of the object. This model was subtracted from the 
original image (i.e. before the first background subtraction). 
This step was done for all selected objects. 

Second background estimation and object mask The out- 
come of the previous step was an image similar to the 
original one, in which large objects have been removed. 
This allowed us to obtain an improved background map, a 
better background subtracted image and a better detection 
of the small objects. In fact, using SExtractor with the 
background-subtracted images after removing the models, 
it was possible to get a careful mask of all objects (ex- 
cept, of course, from those that were removed from the 
image). This new mask would serve to refine the previous 
(manually-done) one. 

Second fit and modeling The fitting process done in the first 
iteration was repeated using a refined mask for each object . 
This helped to get a better fit and, as a result, a better model. 
The main difference with respect to the first iteration is 
that after the subtraction of the model from the galaxies, 
sometimes some residuals remained in the central regions. 
In such cases, these were manually edited using the IRAF 
task imedit. When editing the images only clearly spurious 
residuals were replaced by pixels simulating the background 
signal adding a Gaussian noise whose sigma was computed 
from surrounding regions. This has no effect on an object's 
photometry since the photometry of large galaxies is done 
with another image, avoiding the spurious detection of the 
residuals. 

Construction of the final images In most fields, only two it- 
erations were enough to achieve a satisfactory result. After 
modeling all the selected extended objects and subtracting 
them from the initial image, the background was recom- 
puted again and finally subtracted from the model-subtracted 
image. The output of this procedure was a background- 
subtracted image without the extended modeled objects. 
The final step was to produce the complementary image 
containing only the modeled galaxies. To be consistent in 
the photometry of the extended galaxies we do not use the 
models. We constructed a new image containing the origi- 
nal pixels (background subtracted) of the removed galaxies. 
However, to minimize the effect of the projected objects, the 
intensities of the pixels in which these fell were in fact re- 
placed by the intensities of the models in the same pixels. In 
practice, this was done constructing a new mask of the ob- 
jects contained in the background-subtracted image without 
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Fig.A.2. Comparison of the real counts (continuous line) with 
that coming from the simulations. Long dashes line: Input dis- 
tribution of the simulated stars taken from Besancon mod- 
els. Dotted line: Simulated stars classified as stars. Long-short 
dashed line: Simulated stars detected. 

large objects using SExtractor. Then, the final image was 
computed following these criteria for its pixels: 

- If the pixel did not belong to an object (value in the 
mask equal to zero) then the background-subtracted im- 
age value was kept that corresponded to the original pixel 
of the galaxy. 

- If the pixel indeed belonged to an object (value in the 
mask greater than zero) then the value in the model was 
taken instead of the value in the original image. 

- If the pixel fell in an interchip region then the value from 
the model was used if available. In this way, we were 
able to improve the photometry of large galaxies with 
large regions lost in the interchip regions. 

Of course, for construction, the image of the large galaxies 
also lacks the bright stars. 

This procedure was done in both bands (V,B) removing, of 
course, the same objects. 

Appendix B: Image simulations 

One of the methods most commonly used to check the reliability 
of a procedure of star/galaxy classification is to test the proce- 
dure using images in which synthetic objects have been added. 
Of course, this method relies upon the idea that the artificial ob- 
jects are similar enough (from the point of view of the classify- 
ing program) to the real objects so that one can extrapolate the 
results obtained from the simulations (from which the input and 
the output are known) to the real objects. 

We followed this method using the tasks of the IRAF's pack- 
age artdata. We made catalogs of galaxies and stars to build 
the artificial images that were added to the real images. In this 
way we could test also the effects of the variations of the back- 



ground and other real conditions that are difficult or impossi- 
ble to simulate. The effects of crowding coming from adding 
more objects to real images were measured to produce less than 
10% of the lost objects in the most crowded image. Then we 
proceeded with these images as with the real ones. The result- 
ing parameters from SExtractor were compared between real 
and synthetic objects, showing that they were at first sight quite 
similar, giving us confidence that the extrapolation from simu- 
lations to real objects could be done. However, the comparison 
of the final counts from the real and the simulated objects made 
us distrust the simulations, especially that of the faint stars. The 
reason for this conclusion is illustrated in Fig. IA.2I This figure 
shows the star counts from the original catalogs of WINGS (con- 
tinuous line) and from the simulations. For simulated stars three 
lines have been plotted. The long dashed one represents the in- 
put catalog and the numbers are taken from the models of the 
Galaxy from Besangon. The long-short dashed line represents 
the counts of detected simulated stars. Finally, the dotted line 
shows the counts of the simulated stars classified as stars, which 
should be the counts to be compared with the data from WINGS' 
stars catalog (continuous line). From this figure it can be seen 
that, at the faintest magnitudes, it was easier for SExtractor to 
distinguish a real star than a simulated one, violating the initial 
premise that simulated objects are similar to real ones. 

In the construction of the simulated stars we took bright but 
not saturated stars to have a well sampled point spread function 
(PSF) even at the wings. So, one expects simulated stars to be 
more concentrated than real ones and then easier to be classi- 
fied as stars by SExtractor, which is the opposite of what was 
found. The doubts became greater when simulations were done 
using two different PSFs, one for bright stars computed from 
a bright start and another one for faint stars computed from a 
fainter star. In this case, the fraction of faint simulated stars iden- 
tified as stars by SExtractor increased. Therefore, the simula- 
tions turn out to be too dependent on the input. 

As this seemed to be a problem of the mkobjects task we 
tried with the addstars feature of the daophot package to cre- 
ate the synthetic stars. Since this is a package made to study stars 
we expected a better treatment of the simulations. However, the 
results were similar. 

The origin of the difficulties is not clear. Probably, small vari- 
ations of the local conditions where objects are added produce 
large effects in faint simulated objects. 

All these results convinced us not to use the results of the 
simulations to measure the reliability of our star/galaxy classifi- 
cation at faint magnitudes. However, we did not find such prob- 
lematic behavior at bright magnitudes. Fig. IA. 21 shows that the 
problems of detection of simulated stars (before any classifica- 
tion as galaxy or star) start above V ~ 22 when the detection 
rate drops (see Fig. [5]). For these reason we decided to still rely 
on the simulations to estimate the photometric errors and the de- 
tection rates since, unfortunately, there is no better procedure to 
estimate such quantities. 

Appendix C: Comparison with WFPC2@HST image 

To check the completeness estimations done with simulations 
we performed a comparison using data from the Hubble Space 
Telescope. We downloaded the images of the BCG of A 1795 
taken with the WFPC2 from the HST archiveQ We chose the 
F555W filter which was the one that best matches the V band 

19 Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the data archive at the Space Telescope 
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Fig. C.l. Comparison of the completeness of the WINGS cat- 
alog of the field of A 1795. Upper panel: The continuous line 
shows the completeness computed from simulations while the 
dots with the error bars show the completeness when compar- 
ing from HST data. The error bars are constructing using lcr 
Poisson errors. Lower panel: Number of objects detected in the 
HST image. Since the area is quite small the total numbers are 
also small producing the large uncertainties in the computation 
of the completeness. 



Appendix D: Additional information 

In addition to the data of the single objects found in each field 
we include several tables containing information about observa- 
tional features and peculiarities of either the single clusters or 
the fields on which they are projected. 

Position, redshift, Abell richness, Bautz-Morgan type, X-ray 
luminosity and galactic ex tinction of the whole WINGS sample 
can be found in Table 5 of iPaper j 

In Table ID. 11 we summarize the conversion factors used for 
this work from CCD related units (pixels) to angular unit (arc- 
sees) and from these to linear units (kpc) at the redshift of 
the target cluster using a cosmological model with parameters 
H = 75 km&- 1 Mpc -1 , Q. M = 0.3 and Q. A = 0.7. We have also 
included the effective total area of each image which is the real 
area used to make the catalogs which is slightly smaller than 
the total field of view (for WFC@INT images, the effective area 
is ~ 90% of the total field of view while for WFI@ESO this 
value is ~ 95%). We report the angular sizes of the apertures of 
R = 2Kpc, 5 Kpc and 10 Kpc used to construct our catalogs. 

In Table ID. 21 we list the detection limits in surface bright- 
ness, or surface brightness thresholds, (fiv(Threshold)) as well 
as the V band magnitude at which the detection rate goes down 
to 90%, 75% and 50%. These last values are average values ob- 
tained from the simulations. 

Table ID . 3l reports the positions of the brightest cluster galax- 
ies. In most cases we preferred the coordinates of the peak of the 
emission instead of the coordinates of the barycenter because the 
latter are more affected by irregularities in the outer isophotes. 

Finally, Table lD.4l lists comments or issues about the clusters 
and the fields which we find interesting or useful when working 
with our catalogs. 



images that we used to construct our catalogs. From the mo- 
saic of the WFPC2 we removed the PC chip in which the BCG 
was located since this produced problems for SExtractor and 
we were interested in knowing the completeness at faint mag- 
nitudes. After that, SExtractor was run on the image and the 
resulting catalog was matched and compared with the WINGS 
catalog. We performed the matching against the global WINGS 
catalogs, i.e. including stars, galaxies and objects of unknown 
classification. Fig. IC.ll (upper panel) shows the completeness 
computed for the A 1795 field using the simulations (continu- 
ous line) and the completeness compared with the HST data. We 
also include the errors in the computation of these data since the 
area is quite small and therefore the number of detected objects 
(lower panel) is also quite low. 

Although the HST image is much sharper than our ground- 
based image, it is not much deeper and it also shows problems 
of completeness in the range of comparison (lower panel of 
Fig. IC. 11 1. This and other issues in the matching procedure (such 
as pairs of objects that are not resolved in the WINGS image) 
introduce uncertainties in the comparison so this should be con- 
sidered as a complementary check to the completeness computed 
using the simulations. 



Science Institute. STScI is operated by the Association of Universities 
for Research in Astronomy, Inc. under NASA contract NAS 5-26555. 
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Table D.l. Useful parameters of the WINGS' clusters sample 
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Continued on next page 
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Table D.l. Useful parameters of the WINGS' sample (continued). 
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DM: Distance modulus. 

Total Area: This is the effective area used to construct the catalogs. To compute it the interchip regions have been subtracted. 
DAper '■ Conversion of the physical apertures to angular sizes. 



16 J. Varela et aL: WINGS II: Deep Optical Catalogs 

Table D.2. Completeness limits and surface brightness detection limits. 
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Table D.3. Coordinates of the emission peak of the brightest cluster galaxies 
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Table D.4. Remarks about the individual fields. 



Cluster 



Notes 



A85 

A119 

A133 

A147 

A151 

A160 

A168 

A193 

A311 

A376 



Two very bright stars in (a, <5)v 2 ooo.o=(15.6155°,-21.6094°) and (15.5469°,-21.6080°) 
Very bright star in (16.9983°, 1 .9926°) 
Background structure at z ~ 0.096 
Background group at z ~ 0.060 

Merger of a spiral dominated cluster and a cD dominated one 



Strong spatial variation of the PSF along the declination direction 
Bright star in (41.1317°,36.7629°) 



A500 
A548b 
A602 
A671 



Interacting central pair of galaxies 
Cluster with two similar dominant galaxies 



A754 Two bright stars at (137.1675°,-9.3841°) and (136.9504°,-9.8537°) 

A780 Bright stars at (139.5243°,-12.3193°) and (139.3175°,-12.0079°) 
BCG contaminated by two stars that are at ~ 1 arcmin 
High contamination from a star outside of the field of view. 

A957x Very bright star at (153 .6394°,- 1. 1 178°) 

A970 

A 1069 

A1291 Bright star in (173.0842°,56.0965°) 

A1631a Contaminated by a foreground group at z ~ 0.0144 

The brightest galaxy in the center of the field is a foreground galaxy 
A 1644 Merging between two clusters 

Bright star close to the center (194.3967°,-17.3951°) 



A1736 Cluster without a clear center and with two dominant galaxies quite far away one from the other (> 14' > 700 kpc). 

Possible merger ? 

A 1 795 Bright star at (207 .0113°,26.5154°) 

Bright star at (207.2009°,26.6144°) at only ~ 1.'67 of the BCG 

A1831 

A1983 Two dominant galaxies 

A1991 Bright star at (223.4967°,18.4998°) 

A2107 Only field observed using high-dithering technique 

A2124 

A2149 Very bright star at (240.1033°,53.9671°) 

A2 1 69 Bright spike near the edge of one chip from (4050, 1 ) to (4050, ~ 1 00) 

A2256 Field with a high star density 

Bright star at (255.8912°,78.6300°) 

A2271 Bright star in (257.0958°,78.8285°) 

A2382 Bright star at (328.1596°,-15.8079°) 

A2399 The brightest galaxy in the field is indeed a foreground (z = 0.017) spiral galaxy 

A2415 Bright star at (331.5222°,-5.3578°) 

A2457 Contaminated by a galactic nebula 

A2572a Central pair of galaxies in interaction 

Bright star at (349.4753°,18.7029°) 

A2589 Background group at z ~ 0.17 

A2593 Bright star at (35 1.0991°, 14.4786°) 
A2622 

A2626 Contaminated by A2625 (z = 0.0609) 

A2657 Field contaminated by a galactic nebula 

A2665 Bright stars at (357.6523°,6.0191°) 

A2717 Bright star at (0.8530°,-35.9656°) close to the BCG (~ 2:9) 
A2734 

A3 128 Bright star at (53.0046°,-52.4743°) 
A3158 

A3164 BCGlostin interchip region ! ! 

A3266 Bright star at (67.9736°,-61.3848°) 



A1668 



Continued on next page 
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Cluster Notes 

A3376 Bright star at (90.1511 o ,-39.9517 o ) 

A3395 

A3490 

A3497 About half of the BCG falls in an interchip region 
A3528a Double system together with A3528b 

Member of a group of clusters with A3528a, A3530 and A3532 

Bright stars at (193.5528°,-29.4042°) and (193.4720°,-29.3077°) 
A3528b Only 460 s of total exposure time in V and 180 s in B 

Double system together with A3528a 

Member of a group of clusters with A3528b, A3530 and A3532 
A3530 Member of a group of clusters with A3528a, A3528b and A3532 

BCG in interaction with a close galaxy 

Bright star at (193.8856°,-30.0642°) 
A3532 Member of a group of clusters with A3528a, A3528b and A3530 
A3556 Cluster with two dominant galaxies 

A3558 Bright stars at (202.0540V3 1.5522°) and (201.6474°,-31.4513°) 

A3560 The central galaxy of the field which is also the brightest galaxy in the field is a foreground galaxy (z = 0.0124) 

Bright star at (203.063 1°,-33.2489°) 
A3562 Total exposure times of 180 s in each filter 

Very bad seeing (FWHM t (V) ~ 2"38) 
A3667 Bright star at (303 .24 1 8°,-56. 8468°) 

A3716 Brightest galaxy in the field is a background galaxy (z = 0.0557) 

Two bright central galaxies maybe in interaction 

Bright star at (312.8122°,-52.6243°) 
A3809 Bright star at (326.4170°,-44.1396°) 

A3880 Bright star at (336.9604°,-30.6026°) at ~ 1 !82 from the BCG 
A4059 

IIZW 1 08 BCG in clear interaction 

Very bright star at (318.1888°,2.6429°) 
MKW3s Two bright stars at (230.2168°,7.6827°) and (230.5907°,7.8160°), respectively 
RX 58 BCG with two nuclei separated by ~ 7" 

In the catalog they are considered as two objects 

Two bright stars (14.6187°,26.7859°) and (14.6623°,26.8584°) both closer than 5 arcmin to the BCG 
RX1022 Problem with photometry found, probably extinction by clouds. Corrected using SDSS data. 
RX1740 Field with high density of stars. 

Bright star at (265.3072°,35.7383°) and (265.4106°,35.3667°) 

Z1261 

Z2844 Very bright spike in the border of one chip (150.98 10°,32.6458°) 
Z8338 Field with high density of stars 

Z8852 Very bright star at (347.5279°,7.3640°) 



